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ABSTRACT 


Haystack Pointing System consists of hardware and software which points the 
Haystack 120-foot X-band antenna dish with great accuracies. ‘The Satellite 
Acquisition program, described in this report, generates acquisition scans and 
searches for target returns. Once the target has been sighted, the program tracks 
the target by conical scanning. After acquisition, time correction can be made in 


the orbit computations. 
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HAYSTACK POINTING SYSTEM: SATELLITE ACQUISITION 


INTRODUCTION 


One of the objectives of the Haystack facility is to track satellites for various 
experiments. Monopulse autotrackers have been built for this purpose. The auto- 
trackers require that the antenna beam be pointed at the satellite, before automatic 


> — 


tracking can take place. The satellite program of the Haystack Pointing System 
is capable of computing azimuth and elevation commands from orbit parameters; but 
past experience has shown that the accuracy of the available orbit parameters is too 
poor to place the satellite within the antenna beam. An initial search has to be con- 
ducted to locate the satellite before autotracking can take place. 

An acquisition program has been written for the Univac 490 computer to 
facilitate a more orderly search and acquisition. The primary objective of this pro- 
gram is to find the target and point the antenna at it long enough for the autotracker 
to lock on and follow the target from the acquisition point. However, the program is 
capable of independently tracking the satellite by continuously sensing the discrepancy 
between the computed commands and the actual target position. 

The program was written to be used with the first Haystack radar equipment 
and the sequential doppler processor built by W. F. Relley-. This equipment gener - 
ates a pulse, called a RP2 pulse, whenever the receiver output exceeds a given 
threshold. These pulses can be generated by noise and false targets as well as the 
target which is being tracked. 

Indeed, the program only requires that a RP2 like signal be given the computer 


whenever the target is within the antenna beam. Hence, any equipment that generates 


a A. A. Mathiasen and J. D. Drinan, editors, "Haystack Pointing System: 
Satellite," Lincoln Laboratory Technical Note 1965-36, (9 September 1965). 


va A. A. Mathiasen, editor, "Haystack Pointing System: Mathematical Develop- 
ment for Satellites and Belts," Lincoln Laboratory Technical Note 196549, 
(23 September 1965). 


or. A. A. Mathiasen and J. D. Drinan, "Haystack Pointing System: Control 
Structure," Lincoln Laboratory Technical Note 1966-10,(March 9, 1966). 


4. H. G. Weiss, "The Haystack Experimental Facility, “ Lincoln Laboratory 
Technical Report 365, (15 September 1964). 


such a pulse can be used in conjunction with the acquisition program. R. Silva has 
constructed a CW monopulse autotracker which generates the RP2 pulse whenever 
signals exceed a threshold. The acquisition program has successfully acquired one 
of the Lincoln experimental satellites, LES II, using this equipment. 

The Haystack computer can send pointing commands to the Westford antenna 
over telephone lines’. A. Dockrey has built circuitry at the Westford site to allow a 
RP2 pulse to be sent back to the Haystack U-490 computer under various modes of 
Westford equipment Operation. The acquisition program has acquired satellites by 
means of Westford equipment in the radar mode as well as in the CW mode. 

Since the doppler and range predictions are usually within the required 
accuracies, the acquisition program only searches and acquires in azimuth and ele- 
vation angles. No provision has been made for range or doppler acquisition. 

The program acquires in two phases. First, a search scan is generated which 
grossly looks for the target. It does this by superimposing a long and narrow scan 
along the orbit and searching for RP2 pulses. Secondly, when a hit is received 
(designated by an RP2 pulse), the values of azimuth and elevation for the target are 
observed and a local scan is initiated around these coordinates. The local scan con- 
sists of a set of concentric circles. The radius of the consecutive circles being 
increased by about one beamwidth from one circle to the next. The local scan is 
continuously computed, even after the antenna is controlled by the analog tracker; so 
that if at any time the analog tracker should lose the target, the computer will have a 
correct pointing angle already computed in the buffer region. By simply changing the 
antenna mode to computer pointing, the target will again be acquired. 

When the program has made three successful local scans, a message is printed 
to indicate that the target has been acquired. At this time the operator can ask the 
program to compute and affect a time correction in orbit computations which is 
expected to correct most of the error in predicted angles. If reacquisition is later 
required, the acquisition program will not have to search quite as large an area. A 


local scan may be sufficient for reacquisition after time correction has been made. 


Ds J. E. Gillis, 'Haystack-West Ford Intersite Coupling Link," Lincoln 
Laboratory Group Report 1964-25, (14 May 1964). 


The program is under complete control of the operator. The automatic 
sequence of scans just described can be overriden by means of the keyboard. For 
instance, at any time the search scan can be made to start over, fixed biases can be 
added to the predicted commands, the program computed errors can be set to zero, 
only local scan can be requested if predictions are known to be sufficiently accurate, 


Gic. 
OPERATOR INTERVENTION 


Figure | shows a typical on-line record of operator action. Typical acquisi- 
tion scans are shown by Fig. 2 where the command elevation is plotted as a funetion 
of command azimuth. When the pointing system is requested to point at a satellite, 
the question "Do you want acquisition” is asked and three choices are given as shown 
in Fig. 1. The only other method of getting to the acquisition program is through 
the attention symbol as shown by Fig. 1. When the program is reached through the 


attention symbol, the operator is offered more choices. 
ee Stop Acquisition 


This choice stops the acquisition scans and causes the program to ignore the 


RP2 interrupts. 
va Search Scan 


This choice sets up the normal acquisition mode. The program will first 
scan in a long and narrow area Of uncertainty along the orbit. The scan shown in 
Fig. 3 is superimposed on the computed pointing commands. ‘The length and width 
of the scan is requested through the keyboard. When RP2 pulses are observed, the 
program automatically switches to local scan. When RP2 pulses disappear for a 


sufficiently long period the program goes back to search scan. 
oe Local Scan 


This choice allows the program to superimpose only a local scan on the 
pointing commands. Figure 4 shows a local scan. The RP2 pulses are observed as 


usual. 
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Fig. |. Example of on-line printout. 
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Fig. 3. Search scan. 
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Local scan. 


Fig. 4. 


4, Clear Biases 
This choice clears the computed error biases in the local scan. 
i Set Biases 


This choice sets the azimuth and elevation local scan biases to the keyboard 


entered values. 
6. Time Correction 


This choice causes a time correction to be made in the ephemeris computation, 


based on the errors in cOmmands at that particular time. 
re Haystack Acquisition 


This choice makes the program operate with Haystack RP2 pulses (channel 8 
external interrupts). Seareh or local scan must be requested before the new site RP2 


pulses will be used. 
8. Westford Acquisition 
This choice makes the program operate with Westford RP2 pulses (channel 12 


external interrupts). Search or local scan must be requested before the new site RP2 


pulses will be used. 
PROGRAM INPUTS AND OUTPUTS 


Basically the acquisition program computes four values of azimuth and eleva - 
tion using predicted pointing azimuth and elevation from common storage and puts the 
four resulting values into common storage. However, additional inputs and outputs 


are used. 


Core Storage Inputs 


W(CAZIM):- Predicted satellite azimuth coordinate for i. + 4 seconds in 
revolutions B27. Le i + 2 is the two second time interval for which the next data 


is to be interpolated. 


W(CELEV):- Predicted satellite elevation coordinate for oo +4 seconds in 


revolutions B27. 


W(AZIMOVER):- A code to indicate where the first azimuth point is to be 
placed. When AZIMOVER is positive first azimuth will be placed in the main antenna 


zone. When AZIMOVER is negative first azimuth will be placed in the overlap zone. 


W(TIMEMODE):- A code to indicate if system is running in real time. When 
TIMEMODE is positive, the system is running in real time. Negative TIMEMODE 


indicates Simulated time. 


W(INAZIMADD): - INAZIMADD gives the input azimuth buffer locations. Upper 
half gives the first word address of presently read into buffer. Lower half gives first 
word address of previously read in buffer. 

W(113):- Azimuth in buffer control word. 

Interrupt Inputs 


a 


The program uses external interrupts on channel 8 (range channel) for 


Haystack gear and channel 10 (intersite coupling channel) for Westford gear. 
Keyboard Inputs — 


The program accepts search scan length and width in degrees whenever search 
scan is requested. It also accepts the local scan azimuth and elevation biases when 


operator requests that option. Keyboard entered parameters are in degrees. 


Core Storage Outputs 


W(ACQAZIM), + 1, + 2, + 3:- Four values of azimuth which have acquisition 
commands superimposed on predicted azimuth commands. Values are in revolutions 


bz /. 


W(ACOE LEV), + 1, 4+ 2. +35> Four values of elevation which have acquisition 
commands superimposed on predicted elevation commands. Values are in revolutions 


B27. 
W(TIMECORR): - Time correction in days B28. 


W(TRACKINDIC): - A code which is set to —0O whenever RP2 pulses are 


observed and to + 0 when no RP2 pulses are observed. 


On-Line Printer Output 


When a target has been acquired, a message “Target Acquired" is printed on 


the high-speed printer by means of the printer log program® 
Keyboard Outputs 


Miscellaneous typing out is performed by means of INTERCOM’ as shown by 
Fig. I. 


INTERNAL OPERATION OF PROGRAM 


Initially, the program is made to scan either according to the search scan or 
the local scan. When the target is observed, the equipment external to the U-490 
computer generates a RP2 pulse which appears at the computer as an external 
interrupt. The interrupt answering routine collects the interrupts in a table over the 
two second pointing system cycle. Every two seconds, the acquisition program is 
entered and all of the RP2 pulses are averaged to obtain an average target position. 
Then the predicted position is computed for that precise instant and an error is com- 
puted. This error is added to the future predicted commands and a local scan is 
initiated around these new coordinates. False alarms are allowed in the sense that 
if during the whole local scan no RP2 pulses are received, search will be continued 
by returning to the scan which was previously in operation. If RP2 pulses are observed 
during a local scan, a new local scan will be generated with the new errors as the new 
center point for the scan. Thus, Once the target is observed, the program will track 
the target by continuously making one local scan after another. When three success - 
ful local scans have been made, a message is printed on the high-speed printer and 


only then a request for time correction will be honored. 
Search Scan 


The selection of search scan has received considerable attention. Typical 


prediction errors show that the error along the satellite track is much greater than 


6. J. D. Drinan (editor), ‘Haystack Pointing System: Auxiliary Real-Time 
Programs, " Lincoln Laboratory Technical Note 1966-6, (31 January 1966). 


v A. A. Mathiasen and J. D. Drinan, editors, “Haystack Pointing System: 
INTERCOM, " Lincoln Laboratory Technical Note 1965-39, (9 September 1965). 
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the deviation from the track. These errors specify a long narrow.area of uncer- 
tainty with the highest probability of actual target location occurring in the center and 
along the satellite track. -Obviously, one would like to scan the highest probability 
area first and then proceed to the lower probability areas. Typical constant error 
probability density contours are shown in Fig. 3. 

The search scan also depends on the equipment limitations. The antenna con- 
trol characteristics are such that large pointing errors occur when the antenna is 
asked to follow high rates. This suggests slow scans. Also the radar pulse repetition 
frequency determines minimum cwell time. One would like to have the target in the 
antenna beam long enough to get at least one and preferably several radar returns. 
These considerations dictate a quite slow search procedure. With a constant accel- 
eration scan just slow enough for the antenna servo, an area of one degree by five 
degrees can be scanned in about eight minutes. 

The values of L (the length along trajectory) and W (the width of scan) are 
given to the acquisition program via the keyboard routine. A scan is generated which 
goes along the trajectory L/2 ahead of the predicted value and L/2 behind the 
predicted value, then moves almost One beamwidth (. 04°) in cross-scan direction 
and repeats the scan along the trajectory, then moves almost one beamwidth to the 
other side of the trajectory and repeats again. The lateral distance from the traject - 
ory is increased almost one beamwidth from one along the trajectory scan to the next, 
until a width of W/2 is scanned on each side of the trajectory. If the target has not 
been found then, a new Scan is initiated with parameters L and W increased by 50%. 

It is expected that the satellite acquisition program will be used mainly on 
the horizon when the satellite first comes into view. ‘The first scan simply waits on 
the horizon and adjusts the azimuth so that a satellite could be acquired if it is early 
or late in its orbit. After this first scan, the above described pattern is used. 

The scan along the satellite orbit (S) is computed by using a constant accel- 


eration scan, 
S$ = 2alf2N—|1]]_ (1) 


where a is the constant acceleration, I iS an index which is increased or decreased 
every two seconds, and N is the number of two second intervals in the half scan, and 


is computed by first computing a time (T). 


Aer 


ae (2) 


where L is the scan length and - ee is the maximum allowable acceleration of 


pointing commands. N is then computed such that 
1 & 2Ne=.3 . (3) 


The value of the constant scan acceleration, a, is computed by 
a= ———- ; (4) 


The value of a obtained by this procedure is approximately equal to ia 


From the azimuth and elevation differences performed on the previous values, 


the trackangle, 9, is determined by 
§ = tan |z——,—| | (5) 
] re 


where E. are the elevation values and A, azimuth values corresponding to time 
T., as shown in Fig. 5. pr T+ 2 corresponds to the ume period for which the 
commands are to be interpolated next. 
The elevation command is then computed by adding the elevation component of 
the scan to the elevation predicted position. 
EX = E,+S5 sin 6 
a 6 (6) 
where ES is the elevation coordinate of the command with the scan added. Now a 
azimuth value (A4) is computed so that the resulting point is constrained to fall on 
the satellite trajectory which is indicated on Fig. 5 as the point labeled Xo. This is 


accomplished by a second order fit and extrapolation 
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In order to prevent singularities in the computations, when 9 < 45°, the 
azimuth component of the scan is computed instead by 
AS - A, +S cos @ (8) 


and the elevation is computed to fall on the curve by 


© _E BoMNg =i. nice 


(Eg = Me a2) VE | — Es = hs) 
+ AAR SAT —A_\A es ia: Sees (A ae ©) 


When the elevation is on the horizon, the minimum elevation is used as 


E, and A‘ is computed by Eq. (7). 


The cross-scan is added to the commands by adding the components of the 


cross-scan to the azimuth and elevation commands. 


AC, 


AC eine , (10) 


and 


ra @ AC cos 9 _, (11) 


E 


where AC is the incremental increase in cross-scan from One scan to the next and 
AC, and ACy, are the two corresponding increments, in the azimuth and elevation 
directions respectively. 
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Local Scan 


The local scan configuration is not very critical and does not have to be com- 
puted very accurately, since the scan itself has very small dimensions. The local 
scan which was chosen first points at the center point of the scan for six seconds. Six 
seconds allows the antenna to move into the local scan region and allows all trans - 
lents to be decreased to negligible values. Three circular scans are then performed 
around the center point. Each circle is made in eight seconds, and has a radius which 
increases almost one beamwidth (. 04°) from one circle to the next. The velocity in 
the smallest circle is such that the target will be observed approximately 1.25 
seconds, which allows 25 hits when the pulse rate is 20 a second. The second and 
third circles will allow 12 and 8 hits respectively. Since the scan is relative to 
the predicted target position and moves along in the expected trajectory, the errors 
should be so small that the larger circular scans are never used. This occurs because, 
whenever hits are received during a local scan, the local scan is restarted with the 
scan center at the coordinates where hits were received. 

The circular scans are computed in a novel way. The antenna commands are 
interpolated from four values of azimuth and elevation by means of the interpolation 
program’. The four values given to the interpolation program are such that a smooth 
circle is interpolated between the points. The following values of azimuth and eleva- 
tion are given at the time when new points are to be interpolated between pe and 


TT. #2 
O 


2 
Ist quadrant -R 
2nd quadrant R 
3rd quadrant R 
4th quadrant “Ry 


The value of azimuth and elevation at +. are Py and B, respectively, and R is 


the radius of the scan at the time of computation. 


Se R. Teoste, "Haystack Pointing System: Interpolation," MIT Lincoln Laboratory 
Group Report 1964-57, (28 October 1964). 
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Average Error Computation 


When RP2 interrupts Occur, an interrupt answering routine stores the 
pertinent information in a table, one interrupt at atime. During the two second cycle 
the acquisition program analyses the information in these tables. Each interrupt has 


four words of information as follows: 


bs Azimuth input buffer control word 
as Range command 

Be Doppler command 

4, Azimuth input buffer control word 


In addition, there is a table of weights which indicates the likelihood of the RP2 being 
a true target. Presently the table of weights consists of equal weights, because the 
hardware does not provide measured doppler and range information. 

Every two seconds when the pointing system cycles through the acquisition 
program, the interrupt tables are examined and an average buffer control word 


POW) sag is computed for the previous two second interval by 


1 N 
BCW = au, W. : iZ 
(BCW),,. = x EZ, W\(BCW) (12) 
where N RP2 pulses had been observed and Ww. and (BCW), were the weights and 
the buffer control words of the individual table entries. The average buffer control 
word is rounded to the nearest integer, and the antenna azimuth and elevation angles 
which correspond to this average buffer control word are read. 

Predicted azimuth and elevation angles are also computed for that instant of 
time. Bessel's four point interpolation formula is used for this purpose which is 
similar to the one used in the interpolation program. The azimuth is given by 

Z S 


Z j , 


and 
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Bly) = .%. 2 G..y° EB, 
i=-1 j 


=o ij i (14) 


where A(y) is the value of azimuth interpolated for the time shift y (in fractions of 
2 seconds) as indicated by the average buffer control word. The values of oy are 
a direct result of the Bessel interpolation ecemicisnne. A. are the values of azimuth 
with A, corresponding to the predicted command given at the beginning of the two 
second interval for which A(y) is to be interpolated. Ey) and E. have the same 
meanings for the elevation commands. 

Now the error can be computed which would have occurred, had the predicted 
comimands been given. This is simply obtained by subtracting the interpolated pre- 
dicted commands A(y) and E(y) from the actual values of antenna azimuth and 
elevation as indicated by the average buffer control word and the input angles. 

In the case of the search scan, this average error over two seconds is added 
to the predicted commands and a local scan is generated around these coordinates. 
However, when RP2 pulses are observed during a local scan, the program is allowed 
to finish the circle and then an average error is computed for the whole circular scan. 


This way the antenna corrections are made in the proper direction. 


Time Correction 


If the operator decides that the desired target has been sighted, the program can 
be signaled to make a time correction in the orbit computations. 


The time correction is made by computing two time increments, 


5 6 
A = ee eae am and (15) 
p -] 
Bn Sop (16) 
2 A,-—A ,"* 1r e 
2 -1 
9. D. R. Hartree, "Numerical Analysis, '' Oxford University Press, London, 1955, 


p. 68. 


ai 


whe re E and A, are the antenna pointing angles at the time of time correction 


—A 
p c 
and ED — E. are actually the averaged pointing corrections computed previously for 


and E and a are the predicted target coordinates at the same time. A 


locating the local scan. 


E,—-—E 


Weighting elevation time correction by : = and azimuth time 
A,—A Sy a a 
2-1 : : 
correction by =——.———_,_-,—_;—-. WE get 
Ey = Ej + Ay =o 4 
OE, — E. 5 Ay ~ A.) 
ke es are 7 er rene on 


2 l Z l 


where At is the desired time correction. The azimuth and elevation biases that 


remain will be computed by 


_ At 

AA = Ay (by = 21) (18) 
and 

cK = At 

AE = E,~E,—(£, -E_)) -_ (19) 


Since the new azimuth and elevation corrections AA and AE are approximate, 


it may take a local scan to again precisely point the antenna at the target. 
PROGRAM DETAILS 


The acquisition program has three entries: working entry, initialization entry 
and interrupt entry. The working section is entered every two seconds in the pointing 
system cycle, the initialization section is entered at the beginning of a run or when 
requested through the attention symbol, and the interrupt answering routine is entered 
whenever a RP2 interrupt occurs. Figures 6, 7, and 8 show the working, initializa- 
tion, and interrupt section flow diagrams respectively. The flow diagrams, in con- 


junction with the program listing in the APPENDIX, are self explanatory. 
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Fig. 6. Acquisition working section flow diagram. 
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INTERAZIM 


Update RP2 
Count 


Set up RP2 
Table Location 


Pie. %.. 


3-62-4527 


RP2 
INTERPT 


Save: 
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Set: 
TRACINDIC=-O 
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Update RP2 
Count 


Setup RP2 
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310re: 
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Range, Doppler 


Restore: 
A, Q, B 6 


Interrupt answering routine. 


The initialization section performs a variety of functions. Each function is a 
direct result of the keyboard request, which are explained in the OPERATOR INTER- 
VENTION section. 

The interrupt answering routine serves the function of recording information 
when RP2 interrupts occur. This routine is entered either by the Westford RP2 
interrupts Or the Haystack RP2 interrupts as requested by the operator. The interrupt 
answering routine also sets the TRACKINDIC which is the location of a code to indicate 
that the target is within the antenna beam, and is used by the autotracking program. 

The working section of the program contains most of acquisition program. 
When acquisition has not been requested, most of this program is bypassed. The 
values of azimuth and elevation are corrected for overlap purposes and are passed on 
to the interpolation program. 

The Haystack antenna is capable of travelling 600 degrees in the azimuth 
plane. Figure 9 shows the azimuth travel limits and labels the overlap zones. The 
azimuth angles given the acquisition program are always between O and 360 degrees. 
The acquisition program computes the overlap information so that the antenna goes 
smoothly through north. The detailed computations are outlined by Fig. 10. 

When the acquisition program has been asked to function, after several tests, 
the program computes either a local scan or a search scan. The local scan computa- 
tion is broken into five alternative computations (five values of k), one for each of the 
four quadrants plus an additional one that is used when the local scan is reinitiated 
(k = O case). k=4 is the first quadrant computation; during this time the previous 
circle is examined for RP2 pulses and appropriate computations are made when some 
RP2 pulses have occurred. 

Figure 11 shows the search scan computations logic, while Fig. 12 shows the 
logic when the search scan has to be backed up four points so that continuity is pre- 


served when after a false target the search scan is continued. 
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Pig. 9. Overlap zones. 
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Fig. 10. Overlap computation. 
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Fig. Il. Search scan computation. 
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Fig. 12. Back up search scan. 
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